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!h Abstract 

We utilize classical molecular dynamics simulations to investigate the intrinsic loss mechanisms 
of monolayer graphene nanoresonators undergoing flexural oscillations. We find that spurious edge 
modes of vibration, which arise not due to externally applied stresses but intrinsically due to 
the different vibrational properties of edge atoms, are the dominant intrinsic loss mechanism that 
reduces the Q-factors. We additionally find that while hydrogen passivation of the free edges is 
ineffective in reducing the spurious edge modes, fixing the free edges is critical to removing the 
spurious edge-induced vibrational states. Our atomistic simulations also show that the Q-factor 
•g degrades inversely proportional to temperature; furthermore, we also demonstrate that the intrinsic 
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losses can be reduced significantly across a range of operating temperatures through the application 
of tensile mechanical strain. 



> 

m PACS numbers: 61.46.-w,62.25.+g,68.35.Gy,68.65.La 

m 
oo 
q 

cn 

o 
o\ 
o 



X 



* Electronic address: harold.park@colorado.edu 



1 



I. INTRODUCTION 



Graphene has recently been discovered as the simplest two-dimensional crystal struc- 
ture^. Since then, graphene has been intensely studied for its unusual electrical properties^, 
for its potential in nanocomposites^, and also for its potential as the basic building block of 
future nanoelectromechanical systems (NEMS) 5 | 6 | 7 | 8 | 9 l 

As a building block for NEMS, graphene shows particular promise for ultrasensitive de- 
tection of masses, forces and pressure due to its own extremely low mass. However, the 
key issue limiting the applicability of graphene as a sensing component is its extremely low 
quality (Q)-factor; the Q-factors of a 20-nm thick multilayer graphene sheet were found 
to range from 100 to 1800 as the temperature decreased from 300 K to 50 KP. Recently, 
using multilayer graphene oxide films, higher Q-factors with values up to 4000 were founcP. 
A higher Q-factor is critical to NEMS device performance and reliability as it implies less 
energy dissipation per vibrational cycle, which enables the graphene NEMS to extend its 
operational lifetime by performing near optimal capacity for a longer period of time. Fur- 
thermore, because the mass or force sensing resolution is inversely proportional to Qj 10 * 11 !, 
low Q-factors are the key limiting factors to the development of highly sensitive and reliable 
graphene-based NEMS. 

A recent experimental study by Sanchez et a/P studied the oscillations of suspended 
multilayer graphene sheets in which two edges of the graphene sheet were fixed, while the 
two other free edge surfaces of graphene remained free to oscillate. The Q-factors measured 
in that work were extremely small, between 2 and 30; the low Q's were attributed to extrinsic 
loss mechanisms such as air damping effects. 

Interestingly, the study also found that the free edges of the graphene sheet often had the 
largest vibrational amplitudes during resonance. To verify this, finite element simulations of 
the graphene sheet were performed; by introducing non-uniform stresses in the suspended 
graphene sheet through application of both an in-plane stretch and an in-plane rotation, the 
authors were able to reproduce the large edge modes of vibration observed experimentally. 
However, no correlation between the large edge modes of vibration and the Q-factors were 
established in that work. 

In this letter, we demonstrate using classical molecular dynamics (MD) simulations^ that 
spurious edge modes that arise during oscillation from the free edges of the graphene sheet 
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constitute the key intrinsic loss mechanism in reducing the Q-factors of oscillating graphene 
membranes; extrinsic loss mechanisms 1 ^ such as gas damping and clamping losses are not 
considered in the present work. More significantly, we demonstrate that the spurious edge 
modes arise intrinsically due to the free edges; no externally applied non-uniform stress field 
is necessary to activate these spurious modes of oscillation. We also find that passivation 
of the free edges with hydrogen is ineffective in removing the edge modes, and demonstrate 
that fixing the edge atoms is the most effective way to reduce the intrinsic edge-induced 
dissipation. Finally, we demonstrate the effectiveness of applied tensile mechanical strain in 
dramatically reducing intrinsic energy dissipation across a range of operating temperatures. 

II. RESULTS 

We first show the results of MD simulations that elucidate the fundamental importance 
of edge effects in reducing the Q-factors of graphene nanoresonators. To elucidate the edge 
effects, we studied a 19.7 Ax 127.8 A graphene monolayer that was comprised of 960 carbon 
atoms, and constrained both long edges (127.8 A) to remain fixed during oscillation. We 
utilized three different boundary conditions on the two short (19.7 A) edges of the graphene 
sheet. First, we applied periodic boundary conditions (PBCs), to eliminate any edge effects 
resulting from carbon atoms having a reduced number of bonding neighbors. Second, we let 
the short edges oscillate freely, with no constraints on their motion. Third, we passivated 
the dangling bonds of the short edge carbon atoms with hydrogen. 

We utilized the second generation Brenner potential (REBO-II) 12 for both the carbon- 
carbon and carbon-hydrogen interactions; this potential has been shown to accurately re- 
produce binding energies, force constants and elastic properties of both graphene and hy- 
drocarbons. For all simulations, the graphene sheet was first equilibrated at a specified 
temperature using a Nose-Hoover thermostat!^ for 50000 steps with a 1 femtosecond (fs) 
time step within an NVT ensemble. After the initial thermal equilibration, a sinusoidal 
velocity profile was applied between the two long (fixed) edges of the graphene sheet, where 
the velocity was zero at the long (fixed) edges, and increased to a maximum value in the 
center of the graphene sheet. After applying the sinusoidal velocity profile, the graphene 
monolayer was allowed to freely oscillate for 100000 steps in an energy conserving NVE en- 
semble. The sinusoidal velocity that was applied to induce the oscillations was only 0.036% 
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FIG. 1: (a) Snapshots of the coherent oscillation of a monolayer graphene sheet at 10 K with 
periodic boundary conditions applied, (b) External energy time history showing minimal energy 
loss, (c) Close up of external energy, kinetic energy and total energy. 
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of the total potential energy of the system, to ensure that nonlinear vibrational modes due 
to the applied velocity field would not be present. 

We first show the results at 10 K where PBCs were applied along the short edges of 
the graphene sheet. As observed in Figure (jlji), the graphene sheet oscillates coherently, 
and completely in phase during the entire simulation time. In particular, no difference in 
vibrational amplitude is observed between the edges and the center of the sheet at any 
time. The coherency of oscillation and the lack of dissipation is confirmed by examining the 
external energy time history of the oscillation in Figure (JTJd) , where the external energy is 
defined as the difference between the potential energy during the oscillation and the potential 
energy of the graphene sheet immediately after thermal equilibration and the application of 
the sinusoidal velocity field, but before the actual oscillation of the graphene monolayer has 
begun. Therefore, oscillations of the external energy about the value of eV, as in Figure 
(JTJd) , occur due to the corresponding oscillations of the graphene monolayer that occur due 
to the applied sinusoidal velocity field; any decrease in oscillation amplitude of the external 
energy, as we will show next for the free edge graphene monolayer, results from intrinsic loss 
mechanisms in the graphene monolayer. 

As clearly shown in Figure ([lb), the coherency of the oscillation for the graphene mono- 
layer with PBCs in Figure (Jlji) results in a minimum of intrinsic energy dissipation. Fur- 
thermore, as shown in Figure (lib), the effects of the energy conserving NVE ensemble are 
observed, i.e. the total energy remains constant due to the out of phase oscillations exhibited 
by the kinetic and potential energies. 

We next show the results where the short edges of the graphene sheet were left free, 
or unconstrained; the edges had a zig-zag orientation. This situation corresponds to the 
suspended graphene multilayers that have been tested experimentally by Bunch et alP, and 
Sanchez et al. 7 ; suspension of the multilayer graphene sheet over a trench resulted in two 
edges of the graphene sheet being fixed (the fixing occurs experimentally through van der 
Waal's interactions with the underlying silicon oxide substrate), and two of the edges being 
free. We note again the low Q-factors reported in both experimental studies. 

The results of the MD simulations for the oscillation of the graphene sheet at 10 K where 
the short edges are unrestrained are shown in Figure ([2]). In comparing the results for 
the PBC case in Figure ([l^), it is clear that the oscillation time history for the graphene 
sheet with free edges in Figure (|2^i) is significantly different. We first notice immediately, 
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FIG. 2: (a) Snapshots of the incoherent oscillation of a monolayer graphene sheet at 10 K with 
two free edges showing mode mixing and beating phenomena, (b) External energy time history 
showing significant energy dissipation and beating phenomena, (c) Close up of external energy, 
kinetic energy and total energy. 
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at t = 0, that the free edges of the graphene sheet have a significantly larger vibrational 
amplitude than the interior of the sheet. Furthermore, we notice at t — that the vibrational 
amplitudes of the free (short) edges are opposite in sign. With an increase in simulation 
time, the spurious edge modes begin to propagate from the edges of the graphene sheet 
into the interior, thereby causing different portions of the graphene sheet to oscillate at 
different frequencies; the snapshots in Figure (|2pi) at t = 146/s and t = 980/s clearly 
illustrate how the edges, as well as different sections of the graphene sheet are oscillating 
incoherently and out of phase. We note that during the MD simulations, the edges do not 
always have opposite vibrational amplitudes at t — 0; due to the statistical nature of the 
Nose- Hoover thermostatP', the initial vibrational amplitude at t — at the edges varies with 
each simulation. However, what is constant in all simulations is that the spurious edge modes 
eventually dominate the oscillation of the graphene sheet, similar to that demonstrated for 
a specific case in Figure (|2p,). 

The effects of the incoherent oscillation induced by the free edge effects on the energy 
dissipation are shown in Figure ^p). In comparing the external energy time history for the 
free edge case in Figure ^p) and the PBC case where there are no dangling bonds for the 
edge carbon atoms in Figure ([lb), we immediately see a significant increase in dissipation 
for the free edge case in Figure ^p)- Furthermore, we can detect a beating phenomenon 
for the free edge case in Figure ^p)] the beating phenomenon manifests itself through the 
external energy whereby the energy does not monotonically decrease with increasing time, 
or the number of vibrational cycles. Instead, the energy demonstrates peaks and valleys 
that do not correspond with the vibrational period of the graphene sheet; these peaks and 
valleys occur due to the incoherent mixing and interaction of disparate vibrational periods 
for the graphene sheet with free edges as seen in Figure (j2^i). 

Further insights can be obtained by comparing the kinetic, external and total energies 
in Figure ^p) as compared to the PBC energies in Figure We emphasize again that 

due to the NVE ensemble, the total energy is conserved in the free edge CclSG, clS shown 
in Figure ((2p), just as it is in the PBC case in Figure However, the external energy 

and the kinetic energy are both observed to decrease dramatically in the free edge case in 
Figure (|2p); the total energy remains constant while the components of the total energy are 
decreasing due to the out of phase oscillation of the kinetic and external energies. 

We compare our results for the oscillation of the monolayer graphene sheet with free 
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edges to that of Sanchez et al. 7 . In that work, they also studied the dominance of edge 
oscillations. However, in order to obtain the spurious edge modes, they externally applied a 
non-uniform stress and to the suspended graphene multilayer sheet in the form of both an 
in-plane stretch and an in-plane rotation. Only after applying these states of stress were the 
spurious edge modes detected in their numerical (continuum finite element) simulations. 

However, as shown in the present work, the spurious edge modes occur naturally through 
the MD simulations without application of any non-uniform stresses. Instead, they arise due 
to the fact that the carbon atoms at the edges of the graphene sheet have fewer bonding 
neighbors, and are therefore undercoordinated with respect to the carbon atoms in the 
interior of the graphene sheet. The lack of bonding neighbors means that the stiffness, 
and therefore the vibrational frequency of the edge atoms differs from the atoms within the 
graphene bulk; this difference in vibrational frequency of the edge atoms was illustrated in 
Figure 

We additionally note that both the free edge and PBC calculations were performed at 
10 K, to minimize the effects of energy dissipation that are known to occur in oscillating 
nanoresonators with increasing temperature^. Because the energy dissipation, and thus 
Q-factor degradation occurs for the free edge graphene sheet at low temperature where 
thermal losses are known to be minimized, we have demonstrated that incoherent vibrational 
states resulting from edge effects are the dominant intrinsic loss mechanism for graphene 
nanoresonators. 

Because we have determined that the undercoordination of free edge carbon atoms in the 
graphene sheet is the key source for intrinsic damping effects in graphene nanoresonators, 
we performed simulations in which the edge carbon atoms were passivated with hydrogen, 
to remove the effects of undercoordination. The validity of this idea is further enforced by 
recent works^^, both of which found that hydrogen passivation stabilizes the free edge 
graphene atoms, leaving the graphene edges essentially stress-free. 

The results of the MD simulations of the oscillation of the graphene monolayer with 
hydrogen passivation at 10 K are shown in Figure Q. There are several noteworthy points 
in analyzing the oscillation history in Figure (j3^i). First, we notice that throughout the 
oscillation time history, the edges oscillate nearly in phase with the remainder of the graphene 
sheet. However, the beating phenomena that was also observed for the free edge case in 
Figure pb) is observed again in the oscillation time history and also the external energy 
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FIG. 3: (a) Snapshots of the oscillation of a monolayer graphene sheet at 10 K with its free 
edges passivated with hydrogen showing beating modes, (b) External energy time history showing 
beating phenomena. 



time history in Figure (|3p). The fact that different sections oscillate with different oscillation 
periods is most clearly observed at t — 980/s in Figure ((3^i), where the multiple minima 
and maxima of oscillation amplitude along the graphene sheet are observed. Therefore, 
while hydrogen atoms can effectively energetically stabilize free edge carbon atoms that 
are undercoordinated, hydrogen passivation of dangling carbon bonds at the edges of the 
graphene sheet does not appear to be an optimal solution to mitigate free edge-induced 
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energy dissipation in graphene nanoresonators. 

Finally, we utilize MD simulations to study the intrinsic energy dissipation of monolayer 
graphene sheets where all edge atoms are constrained not to move. In particular, we study 
circular monolayer graphene sheets that are similar geometrically to the circular multilayer 
graphene oxide sheets that have recently been fabricated and testecP, and where extremely 
high Q-factors with values up to 4000 have been found. In addition, the graphene oxide 
multilayers studied by Robinson et a/P were under tensile stress, which has recently proven 
beneficial in enhancing the Q-factors of both metallic^ and semiconducting nanowire J 19 * 20 * 21 ^ 
however, it was not delineated whether the increase in Q found by Robinson et a/P as 
compared to the suspended graphene multilayers tested earlier^ occurred due to the tensile 
stress, or due to the boundary conditions in which all edges of the multilayer graphene oxide 
sheet were fixed. 

We therefore studied the oscillations of a monolayer, circular graphene sheet with di- 
ameter of 42.6 A that was comprised of 547 carbon atoms, in which all edge atoms were 
constrained not to move. We imposed a sinusoidal velocity profile that was a maximum at 
the center of the sheet and decreased to zero at the circular edges to induce the required 
oscillatory motion, where again the kinetic energy of the applied sinusoidal wave was only 
0.03% of the total potential energy of the system to energy linear oscillatory motion. In 
addition, we applied tensile strains of 1,2,3,4 and 5% to the circular graphene sheet by ra- 
dially expanding the sheet to determine the effects of tensile strain on the intrinsic energy 
dissipation of graphene. 

The external energy time history for the circular graphene monolayer oscillating at 100 K 
and at 0% and 1% applied tensile strain is shown in Figure Q. As seen in Figure Q4^l) , when 
no tensile strain is applied, the graphene monolayer does dissipate energy, and the Q-f actor 
is calculated to be about 12000. However, we note that the energy dissipation is clearly 
smaller than the energy dissipation observed for the graphene sheet with free edges that was 
depicted in Figure (J2p) . Furthermore, we note that the rapid energy dissipation in the free 
edge case in Figures ^p) was obtained after 100,000 time steps at a lower temperature (10 
K); the present simulations were run for 3 million time steps in Figure (j4h) at 100 K, where 
significantly less energy dissipation is observed even at elevated temperature and over a 
time scale that is 30 times as long. These results are clearly indicative of the fact that fixing 
the free edges of the graphene sheet, and thereby eliminating free edge vibrations, is the 
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FIG. 4: External energy for circular graphene sheet at 100 K at (a) 0% applied tensile strain, 
(b) 1% applied tensile strain. Dashed lines in both figures are drawn as a guide to the eye in 
quantifying the energy dissipation. 

most effective way to remove spurious edge-induced intrinsic loss mechanisms in graphene 
nanoresonators. 

We show in Figure ^ the variation in Q-factor for the circular graphene sheet as a func- 
tion of temperature. As can be seen, the Q-factor decreases with increasing temperature, 
i.e. from 273000 at 3K to about 2500 at 300 K. We further find that the Q-factor is found 
to be inversely proportional to the temperature, i.e. Q = 1/T. The exponent of 1.0 on 
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FIG. 5: Q-factor as a function of temperature for an unstrained circular graphene sheet. 

the temperature T differs from the exponent of 0.36 found in previous MD simulations^ for 
the Q-factors of fixed/free carbon nanotubes. The most likely explanation for the different 
damping exponent in the case of graphene is because we have removed all edge modes by 
fixing all edge atoms of the circular graphene sheet. Because of this, and because therefore 
the remainder of the carbon atoms in the graphene sheet have a bulk-like bonding environ- 
ment, the Q-factor degrades inversely proportional to the temperature as would be expected 
for a bulk material. In contrast, the fixed/free carbon nanotubes as studied by Jiang et 
al^ had a free end where the undercoordinated carbon atoms were not constrained; the 
non-bulk bonding environment of the free edge atoms of the nanotube are the likely cause 
of the non-bulk temperature damping exponent as compared to a bulk material. 

We also demonstrate in Figures ^p) the external energy time history for the circular 
graphene sheet at 1% tensile strain. As can be observed, there is very minimal energy loss 
when at 1% tensile strain; applying larger amounts of tensile also leads to very small intrinsic 
energy losses during the oscillation of the circular graphene sheet. This phenomenon was also 
observed for other geometries (square graphene sheets), and a wide range of temperatures 
from 3-300 K. This result for graphene is in line with recent theoretical results on metal 
nanowires^, and experimental results where the Q-factors of SiN, Si and SiC nanowires have 
been elevated nearly an order of magnitude through the application of tensile stress^M!^, 
though the mechanisms causing the increase in Q are different between graphene and the 
nanowires. 

Experimentally, it was argued by Verbridge et alW^ that adding tensile stress to SiN 
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nanowires reduced clamping losses by increasing acoustic mismatch between the resonating 
nanowire and the supporting substrate. In the present work, the increase in Q-factor likely 
results from the fact that increased tensile strain in the graphene sheets mitigates thermally- 
driven fluctuations and local variations in vibrational frequency. Therefore, the present 
simulations show that tensile strain, in combination of the removal of dangling bonds for 
carbon atoms that lie on the edges of the graphene sheet, can be utilized to essentially 
remove all intrinsic loss mechanisms in graphene nanoresonators. 

III. CONCLUSIONS 

In conclusion, we have utilized classical molecular dynamics to study the intrinsic loss 
mechanisms, and therefore the reasons underlying the Q-factor degradation in monolayer 
graphene nanoresonators undergoing flexural oscillations. In doing so, we have determined 
that spurious edge modes of vibration, which arise not due to externally applied stresses 
but intrinsically due to the different vibrational properties of edge atoms, are the dominant 
intrinsic loss mechanism that reduces the Q-factors. We further determined that hydrogen 
passivation is ineffective due to the persistance of the beating phenomena that is also ob- 
served for the free edge, or suspended graphene sheets. Finally, we determined that in the 
absence of applied tensile strain, the Q-factor degrades inversely proportional to tempera- 
ture; however, it was determined that the intrinsic loss mechanisms in graphene can nearly 
be eliminated through a combination of applied tensile mechanical strain and boundary 
conditions, i.e. fixing all edge atoms. 
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